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Abstract 

Chain of kinetic equations for non-equilibrium single, double and s-particle 
distribution functions of particles is obtained taking into account nonlin¬ 
ear hydrodynamic fluctuations. Non-equilibrium distribution function of 
non-linear hydrodynamic fluctuations satishes a generalized Fokker-Planck 
equation. The method of non-equilibrium statistical operator by Zubarev 
is applied. A way of calculating of the structural distribution function of 
hydrodynamic collective variables and their hydrodynamic velocities (above 
Gaussian approximation) contained in the generalized Fokker-Planck equa¬ 
tion for the non-equilibrium distribution function of hydrodynamic collective 
variables is proposed. 
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1. Introduction 


The study of nonlinear kinetic and hydrodynamic fluctuations in dense 
gases, liquids and plasma, in turbulence phenomena and dynamics of phase 
transitions, in chemical reactions and self-organizing processes are relevant 
both on kinetic and hydrodynamic levels of description in statistical theory 

of non-equilibrium processes B 3 3 0,3 3 3 3 3 M, IHI, [13 El Q nil, 
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MM 0, [3 • The non-equilibrium states of such systems are far from 

equilibrium. Therefore the study of both the processes establishing the sta¬ 
tionary states with characteristic times of life and the relaxation processes 
to the known equilibrium states, that are described by means of molecular 
hydrodynamics (^, Mi 23. 25 , is of great importance. An important fea¬ 

ture of theoretical modeling of non-equilibrium phenomena in dense gases, 
liquids, dense plasmas (dusty plasmas) is a consistent description of kinetic 
and hydrodynamic processes (isl . 1^ . 27 . 28 . 2^ and taking into account the 
characteristic short and long-range interactions between the particles of the 
systems. In particular, the non-equilibrium gas-liquid phase transition is 
characterized by nonlinear hydrodynamic fluctuations of mass, momentum 
and particle energy, which describe a collective nature of the process and 
dehne the spatial and temporal behavior of the transport coefficients (viscos¬ 
ity, thermal conductivity), time correlation functions and dynamic structure 
factor. At the same time, due to heterogeneity in collective dynamics of 
these fluctuations, liquid drops emerge in the gas phase (in case of transition 
from the gas phase to the liquid phase), or the gas bubbles emerge in the 
liquid phase (in case of transition from the liquid phase to the gas phase), 
formation of which has a kinetic nature described by a redistribution of mo¬ 
mentum and energy, i.e. when a certain group of particles in the system 
receives a signihcant decrease (in the case of drops), or increase (in the case 
of bubbles) of kinetic energy. The particles, that form bubbles or droplets, 
diffuse out of their phases in the liquid or the gas and vice versa. They have 
different values of momentum, energy and pressure in different phases. All 
these features are related to the non-equilibrium unary, binary and s-particle 
distribution functions (which depend on coordinate, momentum and time) 
that satisfy the BBGKY chain of equations. These problems concern the 
consistent description of kinetic and hydrodynamic processes in heterophase 
systems 32, 33|. 


Therefore, the construction of kinetic equations that take into account 


nonlinear hydrodynamic fluctuations |34 l35l. 13611371. l38| is an important prob 


lem in the theory of transport processes in dense gases and liquids. In partic¬ 
ular, this problem arises in the description of low-frequency anomalies in the 
kinetic equations and related "long tail" correlation functions (^, 41 |. 

The main difficulty of the problem is that the kinetics and hydrodynamics 
of these processes are strongly related and should be considered simultane¬ 
ously. Zubarev, Morozov, Omelyan and Tokarchuk 43 . 251 proposed 
the consistent description of kinetic and hydrodynamic processes in dense 
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gases and liquids on the basis of Zubarev non-equilibrium statistical oper¬ 
ator (i^. 45|. In particular, this approach was used to obtain the kinetic 


equation of the revised Enskog theory |43l. |46| for a system of hard spheres 


and kinetic Enskog-Landau equations for one-component system of charged 
hard spheres from the BBGKY chain of equations. 

Zubarev et al obtained the generalized hydrodynamic equations for 
the hydrodynamic variables (densities of the particle, momentum and the 
total energy) connected with the kinetic eq uation for the non-equilibrium 
one-particle distribution function. Later |26l. ^ , these equations were used 
to study time correlation functions and the collective excitation spectrum of 
the weakly non-equilibrium processes in liquids. 

Obviously, the approach proposed by Zubarev et al jd^, 0, 251 and 
Tokarchuk et al j^, can be used to describe both weakly and strongly 
non-equilibrium systems. At the same time, in order to consistently describe 
kinetic processes and nonlinear hydrodynamic fluctuations it is convenient 
to reformulate this theory so that a set of equations for non-equilibrium 
one-particle distribution function and the distribution functional of hydrody¬ 
namic variables, particle number densities as well as momentum and energy 
densities could be obtained. 

In this contribution we will develop an approach for consistent descrip¬ 
tion of kinetic and hydrodynamic processes that are characterized by non¬ 
linear fluctuations and are especially important for the description of non¬ 
equilibrium gas-liquid phase transition. In the second section we will ob¬ 
tain the non-equilibrium statistical operator for non-equilibrium state of the 
system when the parameters of the reduced description are a non-equilib¬ 
rium one-particle distribution function and the distribution function of non¬ 
equilibrium nonlinear hydrodynamic variables. Using this operator we con¬ 
struct the kinetic equations for the non-equilibrium single, double, s-particle 
distribution functions which take into account nonlinear hydrodynamic fluc¬ 
tuations, for which the non-equilibrium distribution function satisfies a gen¬ 
eralized Fokker-Planck equation. In the third section, we will consider one of 
the ways to calculate the structural distribution function of hydrodynamic 
collective variables and their hydrodynamic velocities (in higher than Gaus¬ 
sian approximation), which enter the generalized Fokker-Planck equation for 
the non-equilibrium distribution function of hydrodynamic collective vari¬ 
ables. 
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2. Non-equilibrium distribution function 

For a consistent description of kinetic and hydrodynamic flnctnations 
in a classical one-component flnid it is necessary to select the description 
parameters for one-particle and collective processes. As these parameters 
we choose the non-equilibrinm one-particle distribntion fnnction fi{x]t) = 
{hi{x)y and distribntion fnnction of hydrodynamic variables /(a; t) = (h(a — 
a)y. Here the phase fnnction 


N N 

= - Xj) = ^5(r-rj)5(p-pj) (1) 

j=i j=i 

is the microscopic particle nnmber density. Xj = (r^, p^ ) is the set of phase 
variables (coordinates and momentnms), N is the total nnmber of particles 
in a volnme V. A microscopic phase distribntion of hydrodynamic variables 
is given by 

3 

/(a) = h(a - a) = JJ - a^k), (2) 

m=l k 

where oik = hk, a2k = Jk, hsk = ik are the Fonrier components of the 
densities of particle nnmber, monientnm and energy: 


N 


N 


zkr. 


hk = 

i=i 


Jk = 


—zkrv 


( 3 ) 


i=i 


N 2 N 

= E [III+i E «(I'^«I) 

i=i 




and Omk = (?2'k, Jk,£k) are the corresponding collective variables. $(|rij|) = 
<h(|r; — Fjl) is the pair interaction potential between particles. The average 
valnes (hi(x))* and {6{a — a)y are calculated by means of the non-equilibrium 
A^-particle distribution function g{x^; t), that satisfies the Liouville equation. 
In line with the idea of reduced description of non-equilibrium states this 
function is the functional 


g{x^;t) = p(. ..,/i(x;t),/(a;f),...). 


( 4 ) 


In order to find a non-equilibrium distribution function g{x^] t) we use 
Zubarev’s method 4^, l4a, l47|, in which a general solution of Liouville equa- 
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tion taking into account a projection procedure can be presented in the form; 

= grel{x^-,t) - j Prel{t'))iLNQrel{x^-,t'), ( 5 ) 

J —oo 

where e —)■ +0 after thermodynamic limiting transition. The source selects 
the retarded solutions of Liouville equation with operator zLtv- = 

exp_^_{— dt'{l — Prei{t'))iLi^) is the generalized time evolution operator 
taking into account Kawasaki-Gunton projection Prei{t'). The structure of 
Prei{t') depends on the relevant distribution function grei{x^]t), which in 
method by Zubarev is determined from extremum of the information entropy 
at simultaneous conservation of normalization condition 

J dTNgreiix^;t) = 1, ^ {dxi, .^.^ ,dxN) ^ dx = drdp, (6) 

and the fact that the parameters of the reduced description, fi{x-,t) and 
f{a]t) are hxed. Then relevant distribution function can be written as 

Qrei{x^]t) = exp — ^(t) — J dxj{x]t)ni{x) — J daF(a; f)/(a)|, (7) 

where da is the integration over collective variables: 

da = dnkdjkdek, drik = dRerik dlmrik, £k = dRe^k dlm^k, 

k 

djk djx,]i djy^]idjz^]ii djx,]i d Rejajjk d Imjjj, kj • • • • 

The Massieu-Planck functional <h(f) is determined from the normalization 
condition for the relevant distribution function 

<h(f) = In J dT]\f exp J dx'y{x; t)hi{x) — J daF{a; f)/(a)|. 

The functions 7 ( 0 ;; t) and F(a; t) are the Lagrange multipliers and are deter¬ 
mined from the self-consistent conditions; 

fi{x;t) = {ni{x)y = {ni{x))l^i, /(a;t) = (d(d - a))* = (d(d - a))*,;, ( 8 ) 

where (...)* = J dT m ■ ■ ■ g{x^] t) and (.. = J dT m ■ ■ ■ grei{x^', t). To find 

the explicit form of non-equilibrium distribution function g{x^] t) we exclude 
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the factor F{a]t) in relevant distribution function and thereafter, by means 
of self-consistent conditions dS]), we have 


Qrelix ]t) = Q;^{x ]t) 




W{a;t) 


( 9 ) 


Here 

W{a-t) = j ^ j dV■t)f{a) (10) 

is the structure distribution function of hydrodynamic variables, which could 
be also considered as a Jacobian for transition from /(a) into space of col¬ 
lective variables rj-k, Jk, £k averaged with the "kinetic" relevant distribution 
function 

= exp I - - y dx 7 (a;;f)hi(x)|, (11) 

^fcin(^) _ In j dTj^exp^^— J dx'j{x;t)hi{x)^. 

Here the entropy 

S{t) = -(In Qrelix^-, (12) 

= $(f) + j dx'y{x]t){hi{x)y + j daf{a]t) 

corresponds to the relevant distribution ([9]). In combination with the self- 
consistent conditions ([8]) , it can be considered as entropy of non-equilibrium 
state. In accordance with ([5]), in order to obtain the explicit form of non¬ 
equilibrium distribution function, it is necessary to disclose the action of 
Liouville operator on grei{x^]t) and action of the Kawasaki-Gunton projec¬ 
tion operator, which in our case has the following structure according to 
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(13) 


Prelit)^' = Qrel{x ]t) / dri^Q'+ / dx 


dQrel{x^-,t) 


da 


dgrei{x^;t) 


a(w7ST) 


j j d{hi{x))P 

xi^j dTNni{x)g' - {ni{x)Y j dVNg'^ 

dTNf{a)g' - f{a;t) [ dT^g' 


+ dx da 


dgrei{x^]t) f{a-,t) d\n.W{a]t) 

x(^j dTi^hi{x)g’ - {ni{x)Y j dV^g'^. 

Next, we consider the action of Lionville operator on relevant distribntion 
fnnction (191); 


iLNgrei{x^',t) = - I dx'j{x-,t)ni{x)grei{x^'‘;t) 


N. 


(14) 


iL 


N-, 


/(«;^) 




W{a; t) 

where fii{x) = iL^hYx). Having nsed thereafter the relation 




iL^f^a) — i-h7v/(nk, Jk, ^k) 

^/(a)hk + ^/(a)Jk + ^/(a)4 
UTl\^ yJk 


where 4 = Jk = 4 = *-hAr4, the last expression in (1141) 

can be rewritten in following form; 


iL 


/(«; t) 


N 


lH(a; t) 


QrlUx^-.t) = 
d f{a;t) 


da ly (a; t) 


d fia;t) 

-7777-^ (15) 


‘9nk iy(a; t) 


+ ^k' 


d f{a-,t) 


gL{x^;t). 


"aJkVh(a;f) dekW{a;t). 

Here we introdnced new relevant distribntion fnnction gi{x^, a-,t) with the 
microscopic distribntion of large-scale collective variables 

/(a) 


giix ,a-,t) = gYYl{x ;t) 


iy(a; t)' 


(16) 
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This relevant distribution function is connected with grei{x^',t) by the rela¬ 
tion 

grei{x^;t) = j daf{a-,t)gL{x^,a-,t) (17) 

and is obviously normalized to unity 

j dT^giix^,a;t) = 1. (18) 

Using then the relation fllbl) . the average values with relevant distribution is 
convenient to represent in following form: 

(• • Yg = j da{.. .)i/(a; t), . .)l = j dT^ ■ g^x^, a; t). (19) 


Now in accordance with fllSp and fllbp we can rewrite the action of the Liou- 
ville operator on grei{x^',t) as follows 


iL^greiY^]t) = - I da I dx'y{x;t)hi{x)gL{x^'‘,a-,t)f{a;t) 


N 


( 20 ) 


da W (a; t) 


d d f{a]t) ^ ^ d f{a]t) 
-7T77-;T + 


lU(a; t) 


aJk W(a; t) 


+£k 


d f{a;t) 
dei^ W (a; t) 


gL{x^,a;t). 


Substituting this expression into ([5]), one obtains for non-equilibrium distri¬ 
bution function the following result: 


g{x^;t) = J daf{a;t)gL{x^,a;t) 

+ j da j dx j dt'e^^^'-^Y{t,t')(l - Prel{t')^ 
xni{x)gL{x^, a; t')f{a; t') 

-[ da [ dt'pd'-^Y{t,t')(l - Preiit'Ywia;^)^^- 


( 21 ) 


, 4 d f{a-,t') ^ , d f{a-,t') 
+Jk ■ TT^TTTT-TTT + ^k: 


aJklU(a;f') de-i,W{a;t')i 


d 

drik W (a; t') 
gL{x^,a;t'). 
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and the corresponding generalized transport equations: 


' d p 9 
dt m dr 


/i(x; t)- f dx'-^^{\r - r'|) ■ 


dp dp' 


92{x,x';t) 



d 


i / 1 /N *9 I /(a; t') 

+ (pne{x,Ka]t,t)- 


dsi^) W{a] t') ’ 




d 


d 


d 


E 


A. 


E 


A. 

dev 


j dx' J da' J dt'e''^'^ '''^4’jni.x',\i,a,a']t,t')f{a]t')^{x']t') 
dx' [ da' j dt'e^^'^ ~'''^(l)en{x' ,k,a,a']t,t')f{a]t')'j{x']t') ( 23 ) 

J J — OO 

+E / <'“' /_ <!■ “• “ • *') ■ w,m$) 

[ da' f 

“ J J-oo 9£q 


j da' j dt'ed'-' d ^.A^. (j).^{\^^c^^a,a']t,t') 

k,q 


M^(a; t') 

dec 


^ i /I / /N ^ 1 /(«;^0 

+^4;(k,q.a,a;M)-^}j|^. 

The generalized transport equations fl22|) . fl23|) include the relevant binary 
distribution function of particles g 2 {x,x';t) : 


g 2 {x,x';t) = (6*2(0;, x'))* = (hi(x)hi(x'))* 


( 24 ) 


= / dTN- 2 gqix^;t) = I dag 2 {x,x';a-,t)f{a;t), 


L/^ J. 


where 


5(2 (x,x';a;f) = / drjv_ 2 ^L(x"; a; f) 
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is the binary relevant distribution function of large-scale collective variables. 
The generalized transport kernels (a,/? = {n,j,e}), that describe non- 
Markovian kinetic and hydrodynamic processes, are non-equilibrium corre¬ 


lation functions of generalized fluxes 

(25) 

In{x; t) = (l - P{t))hi{x), (26) 

Jj(k;f) = (l-F(f))jk, (27) 

4(k;f) = (l-F(t))4. (28) 


Here P{t) is the generalized Mori operator related to Kawasaki-Gunton pro¬ 
jection operator Prei{t) by following relation 

Prel{t)a{x)grel{x^]t) = Qrelix^ ] t) P{t)a{x). 

It should be emphasized that in (|25|) the averages are calculated with dis¬ 
tribution function gL{x^,a]t) flT^ . so that the transport kernels are some 
functions of collective variables ak- In equation (|23|) the functions (called 
hydrodynamic velocities) Vn,]i{a;t), Vjk(a;t), represent the fluxes in 

the space of collective variables and are defined as: 

= j dVNn]iQL{x^= (hk)i, 

Vj,k(a;t) = j drArJk^L(a;^,a;f) = (Jk)i, (29) 

Ve,Vi{a]t) = j dr7vdk£)i,(x^,a;f) = (4))^. 

The presented system of transport equations gives consistent description 
of kinetic and hydrodynamic processes of classical fluids which take into 
account long-living fluctuations. 

The system of transport equations ([22]), (|23|) is a not closed due to La¬ 
grange parameter 7 ( 0 :; t), which is determined from the corresponding self- 
consistent conditions. From the kinetic processes standpoint, we must supple¬ 
ment this system of transport equations with the kinetic equation 4 (a;, x'; t), 
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and hence for /^(xi ... x^; f), s < N: 




f2{x,x']t) + iL2/2(x, x'; f) — / dx''{iL{x,x'') + iL{x',x'')}f3{x,x',x'']t) 


= zL2A/2(x, x'; t) — dx"{iL{x, x”) + iL{x', x")}A/3(x, x', x"; t) 


a 



aJv 


, , , , /N ^ f 




/s(xi... x^; f) + iLsfsixi ...Xs] t) 


( 31 ) 


dx.^,^Lix„x.^,)U,{x,...x.,x.^,■t) 



f da f dt'e^^^ l(j)cj{x,x' 

J J —OO ^ 

1/11 i\ 9 ^ /(®; ^0 

+ ^G.£(a;,a; ,k,a;t,f )- 


<9 


) hh(a; t') ’ 

where A/2(x,x';f) = f2{x,x'-,t)-g2{x,x'-,t), Afs{xi ...Xs]t) = /^(xi . ..Xs]t)- 

gs{xi.. .Xs]t). In Eq. (| 3 nl) the two-particle Lionville operator 

iL2 = iLq^x) + iLo{x') + iL{x, x') 
was introdnced. It contains one-particle operator 

iLo(x) = —x = {r,p}, 
m or 
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and also a potential part 


iL{x, x') 




A' 

■ 


Accordingly, in Eg. 0311) . iLg is the s-particle Lionville operator: 


gs{xi ...Xg;t) = {Gg{xi... Xg)y 


dagg{xi... Xg; a; t)f{a; t), 


where 

gg{xi.. .Xg;a;t) = j dTNGg{xi ... a; t) 

is the s-particle relevant distribution function of large-scale variables and 
Gg{xy = ni{xi)....ni{xg). 

Thus we obtained a system of equations for non-equilibrium single, dou¬ 
ble, s-particle distribution functions which take into account nonlinear hy¬ 
drodynamic fluctuations. 

We now discuss the equation fl2^ that is of Fokker-Planck type one 
for non-equilibrium distribution function of collective variables which take 
into account the kinetic processes. The transport kernel in this equation 
(J)nn{,x,x']t,t') describes a dissipation of kinetic processes, while the kernels 
(j)nj{x,k,a]t,t'), (pne{x,k,a]t,t'), (j)jn{x,k,a]t,t’), (j)en{x,k,a]t,t') describe 
a dissipation of correlations between kinetic and hydrodynamic processes. 
To uncover more detailed a structure of transport kernels (f)nn{x] x', a;t,t'), 
(J)Gn{x; x', x", a-,t,t') we consider action of Lionville operator on hi(x) and 
G{x,x') (see Appendix A). That is, taking into account (12^ and (lA.Sp . the 
kinetic equation (l22l) can be written as follows: 
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( 32 ) 


d p 9 

-h — — 

-dt m dr 


X 


fi{x;t)- J dx' j (ia^$(|r - r'l) 
r d d 


dp dp' 


g2ix,x',a;t)f{a;t) = 


dx' f da j dt'e''^'' ■ Djj{x,x'■ -^^{x'-,t')f{a-,t') 


+ / dx' da 


— ■DFj{x,x,a-t,t)- — 


d f f d d ! , d 

' ' ddjpi^x^ X ^a^t^t ^ ' ^^7 ^p " ^ ^ 

d 




X'y{x';t')f{a;t')— '^2 d(^ dt'd'd' k, a; t, t') • 

^ J J —OQ ^ 

, , , /N ^ I /(«;^0 

In the equation fl23|) the quantities 0jj(k, q, a, a'; t, t'), q, a, a'; t, t'), 

0 ej(k, q, a, a'; t, f), 0 £e(k, q, a, a'; t, t') correspond to the dissipative processes 
connected with the correlations between viscous and heat hydrodynamic pro¬ 
cesses. The set of equations fl2^ . (j2^ . flHUl) . (1311) allows for two limiting cases. 
First, if the description of non-equilibrium processes does not take into ac¬ 
count nonlinear hydrodynamic fluctuations, we will obtain generalized kinetic 
equation for the non-equilibrium distribution function of the particles 


■ 9 ^ p d 
dt m dr 


/i(x; t)- f dx'-^^{\r - r'|) ■ 


d d 


-dp dp' 

t 

= [ dx [ dt'e^'" ~'^^(j)nn{x,x';t,t')'y{x';t'). 


g2{x,x';t) (33) 


In this case, the Lagrange parameter 7 ( 0 :; t) is dehned through the function 
fi{x-,t) and, essentially, the equation (l33|) is the closed kinetic equation for 
fi{x-,t) with transport kernel (memory function) ^nn{x,x']t,t'), in which 
the Mori projection operator depends on time. In the case of weakly non¬ 
equilibrium processes kinetic equation (155]) were obtained in jd^, 491 using 
Mori projection operator method. Calculating function ^nn{x, x';t,t') by 
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means of expansions over density, the kinetic equations with linearized Boltz¬ 
mann, Boltzmann-Enskog, Fokker-Planck integrals were obtained. From this 
)oint of view these results can be related with the results of BBGKY theory 
43 |. In our case the BBGKY chain of equations can be obtained assuming 
that A/s ~ 0, namely fs close to Qs- 

Second, if we do not take into account kinetic processes then we will 
obtain generalized (non-Markov) Fokker-Planck equation for non-equilibrium 
distribution function of collective variables, which can be obtained by the 
method of Zwanzig projection operators or by the method of Zubarev non¬ 
equilibrium statistical operator )5l| : 


k 


A. 

dei, 


(34) 





— OO 




ajq w{a-1') 





d 


/^^(k, q, a, a'; t, t') 


de 


d 

dsq W (a; t') 


t 



Now we give two examples of possible applications of the proposed ap¬ 
proach to study the dynamics of simple liquids by taking into account non¬ 
linear hydrodynamic fluctuations. One of these problems is the consistent 
description of kinetic and hydrodynamic fluctuations in turbulence phenom¬ 
ena in gases and liquids. The method of non-equilibrium statistic al op erator 
for the description of turbulent flow was proposed by Zubarev |52 [ l53 | . Main 
idea was to separate the system with the turbulent flow into two subsys¬ 
tems depending on the scale and the character of the motions in the same 
space region. One of these subsystems, that corresponds to large scale and 
regular movement, is described by hydrodynamic equations. The second 
subsystem, that meets to smaller scale and movements with pulsations, is 
described by statistical methods based on Fokker-Planck equation for the 
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distribution functional of the densities of particle number, momentum and 
energy. The interaction between the subsystems leads to an exchange of mo¬ 
mentum and energy between them. The proposed here approach allows to 
describe self-consistently the kinetic and hydrodynamic processes in calcu¬ 
lations of the distribution functional of densities and s-particle distribution 
functions from the equations (125]) . flSU]) . flHT]) . The need for consistent descrip¬ 
tion of kinetic and hydrodynamic nonlinear fluctuations in the investigation 
of turbulence phenomena in gases and liquids were substantiated by Klimon- 


tovich |54 l4|. The other statistical theory in which turbulence is considered 


as non-equilibrium phase transition was proposed by Zubarev, Morozov and 
Troshkin j^. 

Another interesting example might include using the proposed approach 
for consistent description of kinetic and hydrodynamic processes to study the 
phase transitions in non-equilibrium heterophase systems, in particular the 
liquid-gas (vapor) or liquid-glass systems. The heterophase states of physical 


systems were investigated in many papers |^, |^, |^, |^, [18|, |^, |32|, |33 


We know that when heated a water to boiling temperature in it formed va¬ 
por bubbles. On the other hand, above 0°C water is a mixture of water 
and ice, which are present in macroscopic quantities but are not observed 
due to the rapid fluctuations in the system. Namely, the presence of ice 
with structure different from the structure of water (liquid) can explain the 
anomaly density of water in the temperature range from 0°C to 4°C. They 
are the real heterophase systems in which bubble embryoses, drops or small 
crystals have a kinetic nature caused by nonlinear fluctuations, changes in 
temperature, pressure. In our case, the heterophase formations (contain¬ 
ing a hnite particle number in one or another phase) can be described by 
non-equilibrium distribution function /s(x®;f). The kinetic processes within 
heterophase formations are described respectively by the kinetic equations 
fl3T]) . in which right side contains the summands that take into account the 
mutual influence of kinetic and hydrodynamic processes. Obviously, such 
heterophase formations form and decay (with hnite lifetime), by exchang¬ 
ing with particles and energy with the surrounding particles on background 
nonlinear hydrodynamic huctuations of densities of particle number, momen¬ 
tum, energy, the contribution of which increases at phase transformations. 
These nonlinear huctuations are described by the Fokker-Planck equation 
(|23|) . In the process of interaction between kinetic and hydrodynamic huc¬ 
tuations in heterophase systems with some particular change of temperature 
and pressure, the self-organizing particle motions might occur due to spon- 
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taneous symmetry breaking. These particle motions with group velocity 
/s(x®; t) = fs{ri — vt, pi,..., Fg — vt, ps) might lead to an automodel (quasi- 
soliton) propagation of heterophase formations in the system. Such processes 
require separate detailed study due to complex calculative problems of kinetic 
and hydrodynamic transport kernels in transport equations. In connection 
to these processes, we would like to draw attention to the article by Klinion- 
tovich (^, in which to a certain extent a consistent description of kinetics 
and hydrodynamics (taking account diffusion processes) for gas-liquid phase 
transition is realized. In the case in which we do not take into account the 
fluctuations of momentum and energy in our equations, we will arrive for 
a similar Klimontovich equations. These studies require separate considera¬ 
tion, in which the structure function W{a-, t) of collective variables and of the 
hydrodynamic velocities Vn,k{a;t), should be calculated. 

In next section we will perform these calculation. 


3. Calculation of structure function W (a; t) and hydrodynamical 

velocities t) 

In the Kawasaki theory of non-linear fluctuations, the structure func¬ 
tion is approximated by a gaussian dependence on collective variables. In this 
case, as it can be seen, the hydrodynamic velocities nz,k(o-;t), I = are 

the linear function of a. Other approach for calculation of hydrodynamical 
velocities vi^]i{a;t) was proposed on the basis of local thermo dynamics j^. 
The resulting expressions are valid obviously at low frequencies and for small 
values of the wave vector, when the ratios of the local thermodynamics are 
valid. Structure function W{a]t) and hydrodynamical velocities u/ kfa; t ) in 
a case of study of hydrodynamic fluctuations were calculated in us¬ 

ing the method of collective variables (^. The basic idea of this approach 
is that the structure function W{a;t) and hydrodynamic velocities 
can be calculated in approximations higher than Gaussian. Next, we apply 
the method of collective variables (s^ . 62 . 63l for calculating the structure 
function W{a-,t) and hydrodynamic velocities nz,k(o;t). 

First, we calculate the structure function W{a-,t) for collective variables. 
To do this, we use the integral representation for h-functions; 


/(a) = / do; exp 


i,k 


0«,kj 


I = 


( 35 ) 
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Next, using a cumulant expansion for W(a; t) one obtains: 


W(a; t) 




( 36 ) 



yEE 9712^’^" (ki, ks; t)a;zi,kii:^«2,k2} 

h,l2 ki,k2 

X exp I ^Z 7 „(a;;t)|, 

n>3 


where 


az,k = az,k - (a/,k)Lrx> du = W dul^dul^^, a;z,k = t^z'’,k “ *^z!k> ^i,-k = c^z*,k, 

/,k 

Dr,{u;t) =Y1 ^n’"”^"(ki,---,kn;t)u;zi,ki ••••Wz„,k„, (37) 

9^n(ki,...,k„;t) = (azi,ki,-..az„,kJll (38) 

are the non-equilibrium cumulant averages in approximations of the n-order, 
which are calculated using distribution ] t) (ITT]) . We present the struc¬ 

ture function VF(a;t) for further calculations in following form; 

W{a]t)= / do; exp { — ivr E <^Z,kflZ,k (39) 

d Z,k 

^2 

-yEE 9 }t 2 ’^^(ki, k 2 ; f)a;zi,kit^Z2,k2 } 

h,h ki,k2 

x(l + B + iB'^ + iB=> + .^^ + lB"+ .,,), 

where B = X]n >3 -Dn(w; t). If in series of exponent ([39]), namely, e^"^3 Dn{oj-,t) ^ 
one retains only the first term equal to unity, one will obtain the Gaussian 
approximation for W{a]t): 

W^{a;t) = / da;exp{zy ^^a;z,khz,k (40) 

d Z,k 

2 

-yEE 97 X 2 ’*^(ki,k2;f)a;zi,kit^Z2,k2}, 

li,l2 ki,k2 
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where 0 Jt 2 k 2 ; t) are the matrix elements of non-equilibrium correlation 
functions; 



{nn)lin 

(^J)Ln 


im2(ki,k2;t) = 









(41) 


and the non-equilibrium cumulant average 


("kii-k)!! = ("kn-k>L - (*k)L„(ft-k)L„. 


(42) 


For integrating over du in (l39|) we should transform the quadratic form in 
exponential expression into a diagonal form with respect to To this 

end it is necessary to find the eigenvalues of the matrix 0411) by solving the 
equation 

det|in2(ki,k2;t) - ^(ki,k2;t)| = 0, 

£^(ki,k 2 ;t) is the diagonal matrix. Further, obtained eigenvalues F^;(k;t), 
/ = 1,..., 5 of the expression 04O|) are as follows: 

W^{a;t) = / dudetW exp | — (43) 

Z,k 

2 

-yEE £'z(k; t)a;zkWz-k|, 

l k 

where new variables ozk, wzk are connected with the old variables by ratio: 


^nk — 


E 


^/k^/n5 


^/k = 


E 

m=l 


^Zm^mk? 


and uim are matrix elements W = 


Ulu 


1^51 


Wl5 


1^55 


Integrand in 


IS a 


(k;t) 


quadratic function Unk and after integrating over dunk we will obtain follow¬ 
ing structural function in Gaussian approximation VF‘^(a;f): 


1 


W°(a] () = exp { - - ^ B, '(k; ()a,i,a,,_ii} 


( 44 ) 


Z,k 


X 


exp { - - E In 71 det E{k] t )} exp | In det W (k; t )}, 
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or through variables ai]^: 

W^{a;t) = Z{t) exp{ - ^ -Ef(k; t)afka/,-k|, (45) 

^ Z,k 


where 

i' 

Z (t) = exp I - In TT det E{k] t )} exp | In det W (k; t )}. 

k k 

Gaussian approximation, corresponding to Kawasaki theory, is obtained 
as follows. Kinetic processes are ignored in the calculations, then the average 
values are calculated with microcanonical ensemble E = const. During the 
calculation (nJ)^^ = 0, {iJ)% = 0 and the matrix has following form: 


9?t2(ki,k2) = 


transition from the variables h, 


{hn)% 

0 

{he)% 


0 


0 



0 


(ki,k2) 

quadratic form 

in lF'^(a; t) 

, k to h, 

J, h 





{Skn--k)E 

jT—. — ^rik, 
(nkn_k)|; 


where hy_ is the Fourier component of generalized enthalpy. After these trans¬ 
formations Gaussian dependence is obtained as W^{a;t) ~ exp{ —-^l. 

The structure function hF'^(a; t) gives a possibility to calculate (159]) in 
higher approximations over Gaussian moments j^: 


W{a-,t) = W^{a;t)exp{'^{Dn{a-,t))G}, (46) 

n>3 


where one presents {Dn{a]t))G approximately as: 

{Ds{a-,t))G = (D3(a;f))G, 
{D4{a-,t))G = {D4{a-,t))G, 
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{DQ{a]t))G — {DQ{a]t))G - -(.03(0; t))g, 


{Ds{a] t))G = {Dsia-, t))G - {Dsia] t))G{D5{a] t))G - t))^, 


{Dn{a;t))G = 


W^{a-,t) 

1 


X 


/l,... ,/n kl,... jkrx 


W^{a; t). 


..., k„; t) are the renormalized non-equilibrium cumulant aver¬ 
ages of order n for the variables a^k- In expression fH6|) the summands are 
with only even degrees over a since all odd Gaussian moments vanish. 

The method of calculation of the structure function W (a; t) can be applied 
for approximate calculations of hydrodynamic velocities ni,k(a; t). We present 
general formula of velocities consistent with fl2^ in following form: 


vi,kia;t) = J dTNai^kg’^l'l{x^;t)f{a) 

and introduce function Vh(a, A; t): 

W{a,X-,t)= [ p,^:7(a;^;f)/(a), 


so that 


vi,k{a;t) = 


d 


In W (a, A; t) 


^l,k=0 


(47) 


S(-iirAi,k 

We calculate the function Vh(a, A; t) using the preliminary results of the cal 
culation of the structural function W{a;t), and rewrite iy(a, A;t) as: 


Vh(a, A;t)= (ITn du exp | - A;,kark| 

Z,k 


(48) 


X exp I - i7r^a;z,k(az,k - ai,k)|^'re?(a^^G)- 
Z,k 

Now we carry out an averaging in fl48|) with using following cu¬ 

mulant expansion: 


W (a, A; t) = / du exp | — iyr E wz,kaz,k 

Z,k 

[Dn{<^', t) + Dn{X', t) + Dn{u, X] t)] 


( 49 ) 


n>l 
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where 


( ' \ti 

Dn{uj]t)= 

kl,...,k77, 

C„(Ai i) = ^ ^ OT<'I‘-.‘"(ki,.. ^, k„i «)A,„b ^ ^ ^ A,„k., 


D„(a<,Ai() = ^-^ 9J!f''.'»(k,....,k„;() 

v>^n ki,...,kn 


X a;zi,ki • • •n;z„_i,k„_i 

* * ■ ^^71 jkri 5 

with the cumulants of following structure: 


im[)’-’'"(ki,... ,k„;f) = (azi,ki, 

^ \t,C 

• • • ^ln,\<.n/ kin"^ 

(ki,..., k„; t) = (azi,ki, 

^ \t,c 

• • • ^In ,kri / kin ’ 

..., k„; t) = n[{n - j) + {j - n 

+ l)dzi,...,z„_i] 

X (dq,ki) • • • 5 nz^_j.^k„_j 5 • • •) d;„_j+i,k„_j_|_u 

^ \t,C 

• • • 5 ^/rnkn/ kin' 


First, we consider a Gaussian approximation for W{a, A; t), namely we leave 
in the exponent of an integrand only the summands with n = 2 and linear 
over Xi^k. 

W^{a,X;t)= / dwexp |f7r^a;«,k«/,k - *7r^(a,,k)fc'nA,,k (50) 

d «,k l,k 

^2 

“ y X] X] 0 ?t2'’'"(ki, ks; t)a;ii,kit>;«2,k2 

hyh ki,k2 

- ^ 93f®^^’^^(ki,k2;t)a;zi,kiAi2,k2|- 

hyh ki,k2 

Then, transforming this expression in the exponent to diagonal quadratic 
form over variables cuz^k, similarly as for W (a; t), after integrating with respect 
to the new variables wz^k, one obtains; 

f f , jj-2 

W^(a,X;t) = / dwexpj - w ^(az,k)Ln^tk - y y~^(k; t)6z,k5z-k(51) 

d Z,k l,k 

— - In 7rdef^(k; t) + Y^ In detW (k; t) |, 

k k 
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where 


^z,k = 


uJij aj,k + y (k; t) A 


i'.k 


and uij, (k;^) and Ei{k]t) do not dependent on Ai,k- Here the cuniu- 


lant OJt! 


mhf 


(k; t) has the following structure: 


ojt; 




(k; t) — (aj^khj',-k)fcin (%,k)fcm(%',-k)fcm- 


(62) 


Now we calculate the hydrodynamic velocities vi^]^{a-,t) in Gaussian approx¬ 
imation according to the formula fH71) ; 

nz,k(a; t) = ^ ^ In W^{a, A; t) 


d{-i7TXi^i. 


^Lk=0 


(53) 




3,3 


Specifically, we consider the particular case when one can divide the longitu¬ 
dinal and transverse fluctuations for collective variables. That is, we choose 
the direction of the wave vector k along the axis of Oz. Thus, one obtains: 

f TT^ ^ 

W°{a;t)= / £iwexp{s!r^td, 1 , 0 , 4 ,—— ^ 9J;^''‘''"(k:i, k;; t)m,,kitOh.ka 

(54) 


^,k 


^15^2=1 ki,k 2 


“Y ^S’^’*'’^'(ki,k2;fVzi,kiWz2,k2}, 

h ,^ 2=1 ki ,k 2 


where k 2 ; f) are the matrix elements of the non-equilibrium cor¬ 

relation functions of longitudinal fluctuations 

im|(ki,k 2 ;f) = 



(«j")K„ 

if^^Ykin 

(j"^)fcZn 

(j'j")«„ 

(j"^OLn 


(«")«„ 

{^^Ykin 


(55) 


and 


(ki,k2) 

2 (ki, k 2 ; t) are the matrix elements of the non-equilibrium correla¬ 


tion functions of transverse and transverse-longitudinal fluctuations 

®i|’-^(k;,,k; 2 ;t) = 


0 

{f^^xYkin 

Y^^Ykin 

0 


(Jx Jx)fcm 

(Jx 

(Jx^)fc*n 


(Jy Jx)fcm 



0 

(^Jx)fcin 

(^Jy)fc*n 

0 
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( 56 ) 


(ki,k 2 ) 










In this case, the hydrodynamic velocities in the Gaussian approximation are 
as follows: 


t) = 
4k = 
4k (g = 

where 

G„(k;t) = 
Gj(k; t) -- 
^^(kG) 


(^k)Ln + -^1 ^(k;t)(a;iink + a;2iJl|, + i:n3i£k)f^n(k; f), 
(Jk)fcm + ^2 + l^22Jk + '^32^k)^j(k; ^), 

' {^^)kin + -^3 ^(k; ^)(|^13% + Ci;23ji[ + a;33£k)^£(k; t), 


-- a;ii(hkh_k)ll + a;2i(j|[h_k)lt 
= Wi2(hkjik)fcjn + '^22(Jkjik)fcm 
= ^^i3(’^k£-k)fcL + a;23(j4-k)fcL 


+ I^3i(^kh_k)fc[„, 
+ t<;32(£kJ_k)fcin! 
+ W33(^k^-k)4n) 


(57) 


(58) 


and uij are the elements of matrix Vh(k;f). As one can be see, the hydro- 
dynamic velocities (l57)) in the Gaussian approximation for hh'^(a. A; t) are 
the linear functions of collective variables Uk, Jk and e^. It is remarkable 
that if the kinetic processes are not taken into account, then = 1 

(• • •)fcm {■ ■ ■)o is an averaging over a microscopic ensemble fk(a); in this 
case the expressions (1571) for hydrodynamic velocities transform into the re¬ 
sults of previous work in which the nonlinear hydrodynamic fluctuations 
in simple fluids were investigated. The collective variable method 
give a possibility to calculate the hydrodynamic velocities in approximations 
higher than the Gaussian one. In particular, the approximation for the Gaus¬ 
sian, based on (H9|) and hydrodynamic velocities (l57)) will be proportional to 
hi.khi'.k, and transport kernels in the Fokker-Planck equation will be the 
fourth-order correlation functions over the variables o/^k- 

It is important that in Gaussian approximation for kF‘^(k; t) and npj^(a; t), 
the Fokker-Planck equation leads to the transport equations for which 

are similar in structure to the case of the molecular hydrodynamics, averaged 
only over ql{x^, a]t) = Qrd t) ^Fuj^r^■ The proposed approach makes 

possible to go beyond the Gaussian approximation for IF(k;f) and Vi^]^{a]t), 
and hence to do the same in the transport kernels in Foker-Planck equation. 
This allows ns to obtain a nonlinear equation system for (a;,k)* . 

It is noteworthy that kinetic equation (l23l) contains a generalized integral 
of Fokker-Planck type with generalized coefficients of diffusion and particle 
friction in the phase space (r, p, f). This region of changes |r| is limited by val¬ 
ues that correspond to collective nonlinear hydrodynamic processes. 
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This means that in regions of limited the processes are described 

by the generalized coefficients of diffusion and friction, and at small 
they are described by generalized viscosity, thermal conductivity and by 
cross coefficients (j)jj(k,q,a,a';t,t'), cj)j^{k,q,a,a';t,t'), (f)^j{k,q,a,a';t,t'), 
((),e{k,q,a,a';t,t’). Correlations between these regions are described by cross 
kernels (f)nj{x,q,a,a';t,t'), (j)ne{,x,q,a,a';t,t'), (/)£„(k, x', a, a'; f, t'), 0in(k,a;', 
a,a'-,t,t'), that are present both in the kinetic equation and in the Fokker- 
Planck equation. The calculations of these kernels is very important because 
they describe the cross-correlations between kinetic and hydrodynamic pro¬ 
cesses. 


4. Conclusions 


Using the method of Zubarev non-equilibrium statistical operator, we 
have developed an approach 51 | for consistent description of kinetic and 
hydrodynamic processes, that are characterized by non-linear fluctuations. 
We have obtained the non-equilibrium statistical operator of non-equilibrium 
state of the system when the parameters of the reduced description are a 
non-equilibrium one-particle distribution function and the non-equilibrium 
distribution function of the non-linear hydrodynamic variables (densities of 
mass, momentum and energy). By using this operator we constructed a chain 
of kinetic equations (of BBGKY type) for non-equilibrium single, double, 
s-particle distribution functions of particles that take into account the non¬ 
linear hydrodynamic fluctuations. At the same time the non-equilibrium dis¬ 
tribution function of hydrodynamic fluctuations satisfy a generalized Fokker- 
Planck equation. 

We proposed a method to calculate the structural distribution function of 
hydrodynamic collective variables and their hydrodynamic velocities (above 
Gaussian approximation) contained in a generalized Fokker-Planck equation 
for the non-equilibrium distribution function of hydrodynamic collective vari¬ 
ables. In the future studies, we will go beyond the Gaussian approximation 
and carry out approximate calculations of kinetic transport coefficients for a 
specific system of interacting particles. 

The proposed approach of the consistent description of kinetic and hydro- 
dynamic processes can be applied for a description of turbulence phenomena 
in dense gases and complex liquids when a kinetics of a certain component, 
such as behavior of neutrons in the scattering processes, must be take into 
account. It is appealing to apply this approach for description of kinetic and 
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hydrodynamic nonlinear fluctuations in dense dusty plasmas j^, in which 
an important problem is to take into account the non-equilibrium electro¬ 
magnetic processes. We plan to use the outlined approach to investigate the 
non-equilibrium processes in the liquid-vapor phase transition that is charac¬ 
terized by non-equilibrium hydrodynamical fluctuations of masses, momen¬ 
tum and energy of particles. The fluctuations describe the collective nature 
of processes and define the time-spatial behavior of transport coefficients of 
viscosity, thermal conductivity and the process of the emergence of liquid 
drops in the gas phase or the gas bubbles in the liquid phase that has kinetic 
nature. 


Appendix A. 

Now we consider action of Liouville operator on hi{x) and G{x,x')-. 


d \ ^ d ^ 

iLNhi{x) = - -j(r, p) + ^ ■ F(r, p), 

or m op 


(A.l) 


iL^G{x,x') = — 


^ ■ —j(r, p)hi(x') - hi(x)^ ■ —j(r', p') (A. 2 ) 
Or m or' m 

Q ~ d '• 

— ■ F(r, p)fii(a;') + ^1(2^)^ ' F(r', p'), 


dp 


where 


N 


(A.3) 


j(r,p) = 

i=i 

is the microscopic density of momentum vector in coordinate-momentum 
space, 

^ , f) 

(A.4) 


F(r, p) = Y^ - ri)5(p - P, 


3/ 


is the microscopic density of force vector in coordinate-momentum space. 
Taking into account equations flA.ip - (IA.4p . for the kinetic transport ker- 
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nels, we obtain; 


0nn(^7 X ^ tyt ^ 


d , , n d 

— ■D^j{x,x,a-t,t)- — 


dr’ 


(A.5) 


d ^ 9 

-—■DFj{x,x,a-,t,t)- — 


d d d d 

— ■ Djpix, x', ^ ^ ■ Dff{x, x', a; t, t’) ■ — 


where 


dr 


Djj{x,x',a;t,t') = / dTNj{x)T{t,t'){l - P{t'))j{x')pL{x^;t'), 


Dff{x, x', a; t, t') = j dTNF{x)T{t, t'){l - P{t'))F{x')pL{x^] t') 

are the generalized diffusion and the particle friction coefficients in the coordinate- 
momentum space. Moreover, 

/dp/V%(...x';M') = %(r.r';M'). 

J dp J dp' Dff{x, x';t,t') = DFF{^,'i^'','t,'t') 

are the generalized coefficients of diffusion and friction. Similarly, we obtain 
the expression for the transport kernel (j)Gn{x]x',x"-,t,t'): 


(f’Gni^X, X , X , Cl, t, t ) 

d 


d d 

■ Djjn{x, X ,x 

d 


(A.6) 


dr ' ^ ^ 5 ^5 0 ■ Q^n 


d d d d 

■ DFjn{,X, x', x", ^ ^ ■ DFnjix, x', x", d] t, t') ■ — 

d I n ! d d I n I ^ 

■ PjFn{Xi X , X , n, t, t ) ■ '7^ ■ D jnF{x^ X ^ X , Cl, t, t ) • 


dp' dr 


dp" 
d 


d d d 

+ ^ ■ DFFnix, x’, x", ^ ^ ■ DFnF{x, x', x", d; t, t') ■ 
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It is remarkable that expression 



in equation (|23|) with (IA.5j) is the generalized collision integral of Fokker- 
Planck type in the coordinate-momentum space. 
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